TNCT-0601 



Spin polarization and chiral symmetry breaking at 

finite density 



Shinji Maedan0 



Department of Physics, Tokyo National College of Technology, Kunugida-machi, 

Hachioji, Tokyo 193-0997, Japan 



Abstract 

We investigate the possibility of the spin polarization of quark matter at zero tem- 
perature and moderate baryon density. The Nambu-Jona-Lasinio (NJL) model 
including interactions in vector and axial- vector channel (coupling constant G2) 
as well as scalar and pseudoscalar channel (coupling constant G\) is used, and 
the self consistent equations for the spin polarization and chiral condensate are 
solved numerically in the framework of the Hartree approximation. Numerical 
calculations show that in the one flavor model the spin polarization is possible 
at finite density if the ratio G2/G1 is larger than O(l). We find that the inter- 
play between the spin polarization and the chiral symmetry plays an important 
role. If the chemical potential of quarks reaches a certain finite value, the spin 
polarization begins to occur because the quark still has relatively large dynami- 
cal mass, which is generated by spontaneous chiral symmetry breaking. On the 
other hand, when the spin polarization occurs, it operates to lower the dynamical 
quark mass a little. If one increases the chemical potential furthermore, the spin 
polarization becomes weaker and finally disappears because the dynamical quark 
mass is reduced by restoring gradually the chiral symmetry. 
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1 Introduction 



The density of matter in cores of neutron stars is of order O(10 15 g/cm 3 ), and the 
neutron star will give us useful information on how the matter behaves with very high 
density. In neutron stars, it is supposed that nuclear matter or quark matter exists. 
One of the distinctive features of the neutron star is that it has a strong magnetic field 
O(10 12 Gauss). In order to explain the origin of the strong magnetic field in neutron 
stars, the possibility of the spin polarization in nuclear matter has been studied by 
many authors pQ ; however, no definite conclusion has been obtained. Recently, instead 
of nuclear matter, several authors investigate the possibility of the spin polarization in 
quark matter J2j [3], which is expressed by axial vector (^757^) = (ip^ip) with S 
being the spin operator of the quark in the relativistic quantum field theory [3]. This 
article is one of such investigations. We consider the possibility of the spin polarization 
of quark matter at zero temperature and moderate baryon density at which asymptotic 
freedom of quarks does not hold. In such a density region, the interplay between the 
spin polarization and spontaneous chiral symmetry breaking of QCD will be important. 

QCD describes the dynamics of quarks and gluons. For the study of neutron stars 
with high baryon density, one needs to understand the properties of QCD at finite den- 
sity, which have been investigated vigorously by use of effective theories of QCD [5]. 
At low temperature and low density, quarks and gluons are confined and chiral symme- 
try is broken spontaneously, while at low temperature and high density, deconfinement 
occurs and chiral symmetry will be restored, furthermore a new phase called color 
superconductivity may be realized [6J [7] . In Ref . [3] , a coexistent phase of spin polar- 
ization and color superconductivity in high-density QCD is investigated. The authors 
of Ref. [3] use the one-gluon-exchange interaction as an effective interaction between 
quarks, which implies that the density is taken to be high so that one can treat the 
coupling constant a s = g 2 /^ very small. The quark mass is treated as the constant 
parameter not depending on the density, and the numerical calculations are carried out 
for several different values of the quark mass. It is shown that the spin polarization 
remarkably depends on the value of the quark mass. Especially, the spin polarization 
does not occur in the limit of zero quark mass, which is shown analytically by the use 
of the self consistent equation in the framework of the mean field approximation [3]. 

In this paper, we investigate the possibility of the spin polarization of quark matter 
at zero temperature and moderate baryon density, at which spontaneous symmetry 
breaking of chiral symmetry plays an important role. When the quark number density 
increases, the quark mass can not be regarded as a constant, because the chiral symme- 
try is gradually restored and the dynamical quark mass decreases [HI E| • How the spin 
polarization is influenced by the dynamical quark mass which decreases due to restora- 
tion of chiral symmetry? Conversely, how the dynamical quark mass is influenced by 
occurrence of the spin polarization? The point we concentrate is the interplay between 
the spin polarization and chiral symmetry. Although the relation between the chiral 
symmetry breaking phase and the color superconducting phase has been discussed [TU] . 
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we do not consider the color superconductivity in the present paper and focus on the 
relation between chiral symmetry and the spin polarization. 

The method of our approach is as follows. Since the color gauge coupling constant is 
not small at the moderate density, one can not use the one-gluon-exchange interaction as 
an effective interaction between quarks. Instead, effective theories of QCD such as the 
Nambu-Jona-Lasinio (NJL) model [11] are often used at that moderate density region. 
The NJL model realizes spontaneous symmetry breaking of chiral symmetry while it 
does not confine quarks; nevertheless the properties of mesons can be well described 
by the NJL model [9]. The model we use is the NJL type model including interactions 
in vector and axial- vector channel as well as scalar and pseudoscalar channel, which 
model is used to deal with (axial) vector meson modes [HI |13l [9]. 

The paper is organized as follows. The model we use is described in the next section 
and the chemical potential fi of quarks is introduced so as to deal with the system with 
finite density. Using the Hartree approximation (the mean field approximation), we 
can obtain the propagator of the quark in the presence of the mean fields of chiral con- 
densate, quark number density, and axial-vector field. In section 3, the self consistent 
equations for chiral condensate, quark number density, and axial-vector field related 
to the spin polarization are derived respectively in the framework of the mean field 
approximation. These simultaneous self consistent equations are solved numerically in 
section 4. The results of the numerical calculations tell us how the spin polarization 
arises or how the dynamical quark mass changes as the chemical potential of quarks is 
varied. Section 5 is devoted to conclusion. 



2 Formulation and the Hartree approximation 

In this section, the model we use is introduced. Assuming the presence of the mean fields 
of chiral condensate, (ipip), quark number density, ( , 07 t V), and axial- vector, ( / 07 5 7 3 ^), 
at finite density, we calculate the propagator of the quark in the Hartree approximation. 

The following NJL type model [5] with the number of colors N c = 3 and the number 
of flavors Nf = 1 is used, 

C = C + C S + C V , (1) 



£ = ip(0 - 



— — 



Go 



The vector and axial- vector interaction term Ly is needed for treatment of vector meson 
or axial- vector meson in the NJL type model. This lagrangian is chiral invariant if the 
current quark mass m is zero. We make use of this lagrangian in order to investigate 
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the interplay between the spin polarization and chiral symmetry, and regard m and the 
ratio G%jG\ as free parameters. 

The chemical potential \i of quarks is introduced to treat the system with finite 
density, C + fiip^ip. Now, we suppose that the following mean fields exist at finite 
density, 

<^>, $ 7 V>, (^757 3 ^), (2) 

where (V'7573'0) — (V'^sV') i s an expectation value of the z-component of the quark's 
spin S; one can choose the z-component as the direction of the spin polarization with- 
out loss of generality. In the vacuum state (// = 0), chiral symmetry is broken spon- 
taneously, (ipip) 7^ 0, and no spin polarization occurs, (V'757 3 V') — 0- As noted in 
the previous section, color superconductivity is not considered in this paper. Using 
the Hartree approximation (mean field approximation), we obtain the lagrangian with 
bilinear form, 

£mf = - m)i/> + - G^YWYV + G 2 (^ 75 7 3 ^)^757 3 ^ 

= r a (Gl%tf, (3) 

where a,f3 — 1,2,3,4 are the spinor indices, and c,d = 1,2,3 are the color indices. 
The G^ 1 is rewritten more compact form as 



(G^)% = [^-m + G 1 (^) + {^-G 2 ^7V)}7 + G 2 (^757 3 ^)757 3 ] (l) cc ' 



^_M + /i r 7° + f/A757 3 l Jl) cc ', 



a/3 



(4) 



with 



M = m- Gi^ljj), 

fj, t = fi-G 2 (^), (5) 
C/ A = G 2 (^757 3 ^)- 

The chemical potential /x is shifted to /i — G 2 (^7 -?/>) in G^ 1 . In momentum space, 



G~ A l 



i + Hr-f-M +f/A7s7 3 ] (1) 
i>-M +^a757 3 1 (1), 



with 



P M = (po,p) = (9o + A*r,9)- 
The quark propagator G^ (G^G^ 1 = 1) is obtained as (see Appendix A) 

G A = (1) [i>-M + U Albl z 



(6) 
(7) 
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" 1 ' {p 2 - M 2 - XJ\f - AU\ (M 2 + p 2 z ) 

x (pa - M 2 + U 2 A )M + (p 2 - M 2 - U 2 A )t>-2Ul p z ( 7 3 

+2 U APzPfl ( l5 Y ) - ^ (P 2 + M 2 - ^ )( 7s7 3 ) 
+2 U A M (po a 12 - Px a 02 + p y a m 

The propagator G A has four energy poles [3j, p = e n (n — 1, 2, 3, 4), 



ei = e_ = 

e 2 = e + = 
£3 = -ei, 



Pi+ JM 2 +p 2 2 -U A 



p 2 ± + JM 2 +p 2 z + U A 



(9) 



where = p 2 x + p 2 , and ei )2 represents positive energy while e 3j4 represents negative 
energy. Peculiarities of the Fermi surface of e n are described in detail in Ref.[3] 



3 Self consistent equation 

In this section, we derive self consistent equations in the framework of the mean field 
approximation, each of which can be calculated analytically in the zero temperature 
limit. These self consistent equations will be solved numerically in the next section. 
The self consistent equations for (ipip), (^W), and (V'757 3 ?/') are expressed with the 
quark propagator G A , respectively, 



{W) = —iTi{G A }, (10) 

$7V> = -*Tr{G A7 °}, (11) 

^757^) = _, T r{G A7 57 3 }- (12) 

They are rewritten in terms of M, //, and U A that are introduced in Eq.(j5]), 

M = m + iG 1 Tr{G A }, (13) 

Mr = p. + i G 2 Tr{G A7 } = M - G 2 p, (14) 

U A = -tG 2 Tr{G Al5l 3 }, (15) 

where we have introduced the quark number density, p = (ip^ip) = — ?Tr{G J 4 7 }. 



These three equations are simultaneous equations for three unknowns, namely the 
dynamical quark mass M, the "shifted" quark chemical potential p r , and the spin 
polarization U A . 
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3.1 Self consistent equation for M. 



When the system with finite temperature T is concerned, one can employ the imaginary 
time formulation [THE)], according to which the q in the propagator Eq.flSj) is replaced 
by the Matsubara frequency, iuj = i(2j + 1)ttT, (J e Z). Let us begin with the 
calculation of Tr{G/i} in the self consistent equation for M, Eq. (|T3|) . In the imaginary 
time formulation, one needs to perform the Matsubara sum over q Q = iuj 
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Tr{G A } 




T E ^{Ga} 



3=-oo) 

oo 



4 M(p 2 -M 2 + U\ 



(p 2 - M 2 - U\ ) 2 - 4 U\ (M 2 + p 2 z 



(16) 



where Po = 9o + A*r = i^j + /V The result of the sum is (see Appendix B) 

d 3 p 



Tr{G A } = *iV c 



(27T) 



m yw+ji + (-1)^} 



n=l,2,3,4 



M 2 +p 2 



(17) 



where / (e n — /i r ) = l/{e /3 ^ n + 1} is the Fermi-distribution ( (3 = 1/T). In the zero 
temperature limit T — > 0, the self consistent equation for M then becomes 



M — m — GiN c 



d 3 p 
(2^)3 



+ E 



n=l,2 



r„>A/M 2 +p 2 

m{ v /m 2 +p|+ (-i) n ^} 



i c n \i M 2 + p 2 z 



(VA 2 + M 2 



Since the terms coming from the negative energy ( £3,4 < ) contribution in Eq.l iTTj) 
diverge in the zero temperature limit, we regularize these terms by 9 (VA 2 + M 2 — \e n \ ) 
in Eq.(HHD with cutoff A. 



3.2 Self consistent equations for /i r , U a- 



First, the self consistent equation for /i r , Eq. (IT4l) . is calculated. As the Tr{G^} in 
Eq. (jT3l) . TrjC^T } is calculated as follows, 



Tr{G Al °} = iN c 



d 3 p 
(2^)3 



E / ( e « - f^) - 2 

n=l,2,3,4 



The self consistent equation for \i r is then 



/d^p 
— — - ^2 0(fi r -e n ), 

n=l,2 



(19) 



(20) 



in the zero temperature limit T — > 0. 

Secondly, the self consistent equation for C/^, Eq. (IT51) . is considered. Tr{G^7 5 7 3 } 
is evaluated as 

Tr{G A7 57 3 } 



-iN r 



d 3 p 



[U A +(-l)»y/M*+ri\ 

2^ ~ 1 / ( € n ~ 

71=1,2,3,4 6n 



(21) 



and in the limit T — > 0, it becomes [3] 

C? 3 J9 



Tr{G J 47 5 7 } = -*A C 



(2tt) j 



n=l,2 e ™ 



+ E 

n=3,4 



{f/ A + (-l)yM 2 +p2} 



f22) 



where the second term in the right-hand side represents the contribution of the Dirac 
sea ( e 3j4 < ). Here, we define the spin polarization of the vacuum state fi = ( 
H r = ) to be zero by subtracting as follows, 



U A = -iG 2 Tr{GU 757 3 } - Tr{G A7 57 3 }L=o 



(23) 



the definition implies that the contribution of the Dirac sea to the spin polarization is 
neglected. The self consistent equation for Ua eventually takes the form 



-G 2 N C 



(2tt) 3 



(24) 



71=1,2 



in the zero temperature limit. 

We finally obtained the self consistent equations for the dynamical quark mass M, 
Eq. (jl8p . the "shifted" quark chemical potential fi r , Eq. (j20p . and the spin polarization 
Ua, Eq. ()24p . in the zero temperature limit. Each of three equations can be calculated 
analytically if < Ua < M, and the results are shown in Appendix C. 



4 Numerical calculations 

In this chapter, the numerical solutions for the three self consistent equations Eq. (fl~8l) . 
(120]) . and (j24p are obtained, by which we can know how the quark spin polarization Ua 
and the dynamical quark mass M behaves as the quark chemical potential \x varies. 
Giving the five input parameters ( Gi, G%, m, A, fj, ), we solve these simultaneous 
equations of three unknowns ( M, /i r , Ua)- 



6 



700 
600 
500 
400 
300 
200 
100 



300 400 500 600 700 
chemical potential 



800 900 



Figure 1: Dynamical quark mass M(MeV) as a function of the chemical potential 
\x (MeV) when the ratio is GijG\ = 0. The current quark mass m is taken to be MeV 
(solid), 5 MeV (short-dashed), and 130 MeV (long-dashed), respectively. 

At the beginning, let us discuss the input parameters. In the NJL model analysis, 
the ratio GijG\ is taken to be G%jG\ ~ 0.7 in the vacuum state when the number of 
quark flavor is threejl3]. We shall, however, regard GijG\ and m as free parameters in 
our one flavor model since we are interested in studying the relationship between spin 
polarization and spontaneous symmetry breaking of chiral symmetry at finite density 
by the use of this effective theory. In other words, the pion decay constant or the vector 
meson mass calculated by our model with these input parameters does not necessarily 
fit the experimental data. Now, before going into the case of G2/G1 7^ 0, let us consider 
the case of GijG\ = 0, which enables us to concentrate on the issue of chiral symmetry. 



4.1 Numerical results with G2/G1 = 

When G 2 = 0, input parameters are ( G\, m, A, \i ), and an unknown is (M) because 
of Ua = and fi r = fi. The value of the effective quark mass (dynamical mass) M 
in the vacuum state \i — is determined by the current quark mass m, the coupling 
constant G\, and the cutoff A. We put G\ and A as follows, 

Gi = 0.95 x 10~ 5 MeV" 2 , A = 900 MeV. (25) 

For several choices of the current mass m = 0, 5 MeV, 130 MeV, the self consistent 
equation Eq.fllEU gives the value M = 295 MeV, 316 MeV, 601 MeV, respectively. In 
Fig.l, the numerical results of the M's self consistent equation Eq.flTSJ) are shown for 
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m — 0, 5MeV, 130 MeV, respectively. These graphs show that, if the current mass 
m becomes heavier, the value of the chemical potential p at which chiral symme- 
try is restored becomes larger. The reason of such a behavior is as follows. From 
the self consistent equation Eq. ffT8l) (see also Eq. flC.ll) ). the dynamical mass M takes 
the constant value Mq when p r is in the range < p r < Mq. If p(= p r ) ex- 
ceeds the value M , the dynamical mass M decreases and the quark number density 

3 

p = N c (p 2 — M 2 ) 2 6(p — M)/3ir 2 begins to take a positive value. 
4.2 Numerical results with G2/G1 > 

All the numerical calculations are carried out with the input parameters of G\ and A 
given in Eq. ([2"o"j) . As mentioned before, m and G2/G1 are regarded as free parameters. 
At first, we put m = and find the value of G2/G1 by numerical calculations so that 
the spin polarization occurs (Ua > 0) in some chemical potential region (p < A). The 
same procedures are performed for other values, m = 5 MeV, 130 MeV, respectively. 

First, we start with the chiral limit case, m = 0. By numerical calculations for 
the three self consistent equations Eq. ffT8l) . Eq. f[20|) . and Eq. (l24"l) . we find that it is 
necessary to satisfy G2/G1 ~ 7.37 if m = in order to occur the spin polarization. 
When G2/G1 = 7.37, the spin polarization Ua has the nontrivial solution in a chemical 
potential range of about 370 MeV to 470 MeV, and the numerical results are shown 
in Fig. 2 and Fig.3. As seen from Fig. 2, chiral symmetry is restored completely at 
p ~ 630 MeV, which value is considerably larger than the value p ~ 340 MeV at 
which chiral symmetry is restored completely in the case of G2/G1 = and m = 
in Fig.l. This is because the lagrangian Eq.([T]) contains the vector type interaction 
and thereby the chemical potential is shifted by an amount G 2 p, p = p r + G 2 p [51 [TB]. 
This graph Fig. 2 has the peculiarity; the line of the graph is lowered in the range 
370 MeV ~ p ~ 470 MeV. This peculiarity will be discussed later. In Fig.3, the 
nontrivial solution of Ua is figured, and one finds that the quark spin polarization 
occurs in the chemical potential range 370 MeV ~ p ~ 470 MeV. It is notable 
that the spin polarization disappears when the chemical potential becomes too large, 
p ~ 470 MeV, and this point will be also discussed in the later subsection. 

Next, when m = 5 MeV, we find by numerical calculations that it is necessary to 
satisfy G2/G1 ~ 6.3 in order to occur the spin polarization Ua > 0. One can see that 
the necessary value G2/G1 ~ 6.3 in the case of m = 5 MeV is smaller than the value 
G2/G1 ~ 7.37 in the chiral limit case m = 0. Although our NJL model of one flavor at 
finite density differs from the NJL model of three flavor in the vacuum state, the above 
value G2/G1 ~ 6.3 is fairly larger than the value G2/G1 ~ 0.7 which value is expected 
in the NJL model of three flavor in the vacuum state [13]. We then take another large 
value of m and find the value G2/G1 necessary for obtaining spin polarization Ua > 0. 

Finally, when the current mass is taken to be m = 130 MeV, we find that it is 
necessary to satisfy G 2 /Gi ~ 1.89 so as to occur the spin polarization Ua > 0. 
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Figure 2: Dynamical quark mass M(MeV) as a function of the chemical potential 
fjL (MeV) when the ratio is G2/G1 = 7.37 and the current quark mass is m = 0. The 
M in Fig. 2 and the nontrivial solution Ua in Fig. 3 are the solution of the simultaneous 
equations Eq.JEJ, Eq.fl2DJ, and Eq.(|24|). 




Figure 3: Spin polarization Ua (MeV) as a function of the chemical potential \i (MeV) 
when the ratio is GijG\ = 7.37 and the current quark mass is m = 0. The M in Fig. 2 
and the nontrivial solution Ua in Fig. 3 are the solution of the simultaneous equations 
Eq.(HBD, Eq.([20D, and Eq.(f24D. 
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Figure 4: Dynamical quark mass M(MeV) as a function of the chemical potential 
fi (MeV) when the ratio is GijG\ = 1.89 and the current quark mass is m = 130 MeV. 
The M in Fig. 4 and the nontrivial solution Ua in Fig. 5 are the solution of the simulta- 
neous equations Eq.flHD, Eq.(|2QJ), and Eq.(|24l). 




Figure 5: Spin polarization [/^(MeV) as a function of the chemical potential fi (MeV) 
when the ratio is GijG\ = 1.89 and the current quark mass is m = 130 MeV. The M 
in Fig. 4 and the nontrivial solution Ua in Fig. 5 are the solution of the simultaneous 
equations Eq.JEJ, Eq.(f20D, and Eq.flM]). 
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While GijG\ ~ 1.89 is still larger than the value 0.7, we understand that the necessary 
value G2/G1 for obtaining spin polarization becomes smaller if the current quark mass 
m becomes larger. The numerical results for m = 130 MeV and G 2 /Gi = 1.89 are 
represented in Fig. 4 and Fig. 5. The value of fi at which the dynamical mass M begins 
to decrease in Fig.4 is larger than that value of [i in Fig. 2, because the current mass m 
has larger value 130 MeV in Fig.4 than in Fig.2. When the current mass m is 130 MeV 
(constituent mass M = 601 MeV), the spin polarization occurs in the chemical poten- 
tial range 730 MeV ~ [i ~ 860 MeV as seen from Fig. 5. When, on the other hand, 
m is (constituent mass M = 295 MeV), the spin polarization occurs in the range 
370 MeV ~ ~ 470 MeV as seen from Fig. 3. The difference between these two ranges 
of fi in Fig. 5 and Fig.3 comes from the fact that the chemical potential fi must exceed 
the value of the constituent mass M so that the quark number density p has a positive 
value. 

4.3 Effect of the spin polarization Ua on the dynamical quark 
mass M, and vice versa. 

We have solved numerically the simultaneous self consistent equations for M, /i r , and 
Ua when the current mass m = 0, 5 MeV, 130 MeV, respectively. Since these equations 
are the simultaneous ones, Ua has an effect on M, and vice versa. In this subsection, 
we study the effect of Ua on M and the effect of M on Ua through the case of m = 0. 
The behavior of the dynamical mass M in Fig.2 is influenced by the spin polarization 
Ua, and the behavior of Ua in Fig.3 is influenced by the dynamical mass M. 

First, let us consider the effect of the spin polarization Ua on the dynamical mass 
M. As mentioned before, the behavior of M in Fig.2 has the peculiarity that the line 
of the graph is lowered in the range 370 MeV ~ p ~ 470 MeV. This peculiarity will 
be caused by the influence of the spin polarization Ua being non-zero in that range of 
/i as seen in Fig.3. In order to confirm this, we solve numerically the self consistent 
equation for M and that for /i r with setting Ua = in the whole range of fj,, < /i < A. 
The numerical result is shown with the dotted line in Fig. 6, and the graph of Fig.2 is 
also shown in Fig. 6 with the solid line for comparison. When the solid line is compared 
with the dotted line in Fig.6, we find that if the quark spin polarization occurs Ua > 0, 
it operates to lower the dynamical quark mass M which is generated by spontaneous 
symmetry breakdown of chiral symmetry. 

Next, let us consider the effect of the dynamical mass M on the spin polarization 
Ua- As seen in Fig.2, M is a monotone decreasing function of fi, which passes through a 
point (n, M) = (452 MeV, 219 MeV). As mentioned before, the behavior of Ua in Fig.3 
has the peculiarity that the spin polarization Ua disappears when the chemical potential 
becomes too large, /i ~ 470 MeV. This peculiarity may be caused by the fact that the 
dynamical mass M is the monotone decreasing function of fi. In order to confirm this, 
we examine how the Ua behaves when the M is taken to be a constant. Specifically 
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Figure 6: Dynamical quark mass M(MeV) as a function of the chemical potential 
\i (MeV) when the ratio is G 2 /Gi = 7.37 and the current quark mass is m = 0. The 
dotted line represents the dynamical mass when there is no spin polarization ( Ua — 
). The solid line represents the dynamical mass with the nontrivial solution Ua (the 
same result shown in Fig. 2). 




Figure 7: Spin polarization L^i(MeV) as a function of the chemical potential \x (MeV) 
when the ratio is G 2 /G 1 = 7.37. The solid line is a part (/i > 452 MeV) of the non- 
trivial solution shown in Fig. 3. The dotted line is obtained when we set the dynamical 
mass M = constant = 219 MeV. 
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speaking, we solve numerically the self consistent equation for Ua and that for \i r with 
setting M = constant = 219 MeV in the range of /x > 452 MeV. The numerical 
result is shown with the dotted line in Fig. 7, which shows that the spin polarization 
Ua becomes larger as \x increases even in the range of large /i ~ 470 MeV if the M is 
constant ( 219 MeV). In Fig. 7, a part of the graph of Fig. 3 ( fi > 452 MeV ) is also 
shown with the solid line. By comparing this solid line with the dotted line in Fig. 7, 
we find that, when the chiral symmetry is gradually restored as the chemical potential 
H increases ( /i > 452 MeV ), the dynamical quark mass M reduces (M < 219 MeV, 
see Fig. 2) and this reduction of the dynamical quark mass M operates to suppress the 
spin polarization Ua (see the solid line in Fig. 7). 

So far we have argued the special case of m = 0, yet the qualitative properties 
of the effect of Ua on M and the effect of M on Ua are the same with the cases of 
m = 5 MeV and 130 MeV. 

5 Conclusion 

We studied the relation between the quark spin polarization and spontaneous chiral 
symmetry breaking at finite density by use of the NJL type model with the number of 
colors N c = 3 and the number of flavors Nf = 1, which includes interactions in vector 
and axial- vector channel. The self consistent equations for chiral condensate (ipip), 
spin polarization (ip'T,^), and quark number density (ip^ip) are derived respectively 
in the framework of the Hartree approximation (mean field approximation). These 
simultaneous equations are solved numerically with regarding the quark current mass 
m and the coupling constant ratio G 2 /Gi as free parameters. Consequently, we find 
that in the flavor Nf = 1 model the spin polarization is possible at finite density if the 
ratio GijG\ is larger than 0(1). When the current quark mass m becomes heavier, 
the necessary value G<ijG\ for obtaining the spin polarization becomes smaller. For 
example, if m = 130 MeV, the necessary value is GijG\ ~ 1.89. 

The spin polarization occurs in the following way. The spin polarization begins 
to occur at a certain value of the chemical potential, at which the dynamical quark 
mass is still a few hundred MeV. The reason why the spin polarization occurs are 
that the quark number density is finite and that the dynamical quark mass remains 
relatively large value about a few hundred MeV. If the effective mass M of the quark 
is taken to be a small value such as 5 MeV, the spin polarization does not take place. 
When one increases the chemical potential furthermore, the spin polarization becomes 
weaker and finally disappears, where the chiral symmetry is gradually restored and the 
dynamical mass M is reduced from its value in the vacuum state. The reason why 
the spin polarization disappears at the large chemical potential is that the dynamical 
quark mass is reduced by restoring gradually the chiral symmetry. The effects of the 
dynamical mass M on the spin polarization Ua are discussed above; conversely, what is 

2 The M — constant = 219 MeV is not a solution of the self consistent equation for M. 



13 



the effect of the spin polarization on the dynamical mass? The numerical calculations 
show that, if the spin polarization occurs Ua > 0, it operates to lower the dynamical 
quark mass M which is generated by spontaneous chiral symmetry breaking. 

In studying the quark spin polarization with the moderate density at which chiral 
symmetry is partially broken spontaneously, it is necessary to consider the relation be- 
tween the spin polarization and spontaneous chiral symmetry breaking [3]; the present 
paper attempts to investigate such issue. However, there remains several problems not 
discussed here, and we show below; 

• It is necessary to make sure that the nontrivial solution Ua of the self consistent 
equations gives the minimum value of the thermodynamical potential. 

• The effect of quark confinement is not considered, since the NJL model we use is 
an effective theory that does not confine quarks. We therefore treat the system 
at finite density described by the NJL model as quark matter in the present 
paper. Although the lattice QCD simulation demonstrates quark confinement, it 
is difficult to apply it to the system with finite density. The development of this 
field is anticipated. 

• We need to study the case of £77(2) or SU(3) flavors in which the condition 
of charge neutrality is imposed. According to the NJL model analysis in the 
vacuum state, the ratio G2/G1 is taken to be G2/G1 ~ 0.7 when the number 
of flavor is three. In our one flavor model at finite density, it is necessary to 
satisfy G2/G1 ~ 1.89 in order to occur spin polarization if the current quark 
mass is m = 130 MeV. We therefore need to study the spin polarization around 
G2/G1 ~ 0(1) in the case of three flavors. 

• If the color superconductivity is considered, one needs to investigate the relation- 
ship among the spin polarization, spontaneous chiral symmetry breaking, and the 
the color superconductivity. 

• We have done the numerical calculations with T = 0. It would be interesting to 
study numerically the finite temperature case T > 0. 
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Appendix 



A The calculation of the propagator Ga> 

The Ga-i Eq-dS]), is calculated as follows. Following the notation used in the textbook 
of Bjorken and Drell [T?| . we write 



By the way, the inverse of the matrix 



can be obtained by [18J 



where 



x =(cd)- < a - 2 > 



U = {B- AC^D)- 1 , 

Y = (C - DB- l A)~\ (A.4) 
provided and C _1 exist. Using this formula, we get Eq.®. 

B Summation over the Matsubara frequency 

Let us use the standard technique [HI [T5] to calculate the sum YfjL-ao over the Mat- 
subara frequency u>j = (2j + l)nT. Let g(p°) be a function of p° which decreases faster 
than l/p° for \p \ — > oo. Putting p = ity + /V, we obtain the following equation by 
the contour integration [15], 



1 v 2, ( ■ _l \ \ " Res 9{Po) 



Now, in Eq. fTBjl . we take as follows 

4 M(p 2 — M 2 + U\ 



g(po 



(p2 _ M 2 - f/2 )2 _ 4 f/2 ( M 2 + p 2) 

4M(p 2 -M 2 + £/*) ^ (-l)" 1 (R2) 



2(ei-e 2 ) ^ e„ (p -e n ) 
then Eq. ffTTl) is derived by the formula Eq. (IB. lj) 
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C Explicit expressions for the self consistent equa- 
tions 



The self consistent equations Eg. (ITS]) . Eg. (l20l) . and Eq. fpM]) with T = can be calcu- 
lated analytically when < Ua < M. The results are 



M = m 



(27T) 



-M 



|(// r - Wa)^{^-Ua) 2 -M* 



+ (2U A ^r-U A )-M 2 )\og 



+|(// r + 3U A )^(^ + U A ) 2 -M 2 



{lh ~ U A ) + /MaH^I 



M 



je(^-(M + u A )) 



+ (-2U A (fir + U A )-M 2 )\og 



(fir + U A ) + V(/ir + U 

M 



M 2 i 



(2tt)* 



■\e(yL T -{M-U A )) 
(C.l) 



fit = /i 

-G 2 



(27T) 



-2M 2 + (> r - ?7 A )(2/. r + Ua)) - U A ) 2 - M 2 



N c 
(2nf 



-U A M Z log 



2i (/i r - C/ A ) + ^ - U A ) 2 - M 2 



M 



6{^-{M + U A )) 



(C.2) 



(27T) 



| 1 (4Af 2 - - U A )(^ + 2 E^)) ^ - U A ) 2 - M 2 



, , rMwll (/i r -^) + 1 /(/i r -^) 2 -M 2 1 
-(/i r - 2 f/ A )M 2 log 



.1/ 



-G, 



(2vr) 2 



|0(^-(M + C^)) 
(C.3) 



16 



References 

[1] I. Vidafia, A. Polls, and A. Ramos, Phys. Rev. C65, 035804 (2002); I. Vidafia and 
I. Bombaci, Phys. Rev. C66, 045801 (2002), and references therein. 

[2] T. Tatsumi, Phys. Lett. B489, 280 (2000); A. Niegawa, Prog. Theor. Phys. 113, 
581 (2005). 

[3] E. Nakano, T. Maruyama, and T. Tatsumi, Phys. Rev. D68, 105001 (2003); T. 
Tatsumi, T. Maruyama, and E. Nakano, Prog. Theor. Phys. Suppl. 153, 190 (2004). 

[4] T. Maruyama and T. Tatsumi, Nucl. Phys. A693, 710 (2001). 

[5] For a review, see M. Buballa, Phys. Rep. 407, 205 (2005). 

[6] B.C. Barrois, Nucl. Phys. B129, 390 (1977). 

[7] D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 

M. Asakawa and K. Yazaki, Nucl. Phys. A504, 668 (1989). 

[9] For a review, see S. P. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 



[10] For example, M. Huang, |hep-ph/ 0409167, 

[11] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 

[12] M. Takizawa, K. Tsushima, Y. Kohyama, and K. Kubodera, Nucl. Phys. A507, 
611 (1990). 

[13] S. Klimt, M. Lutz, U. Vogl, and W. Weise, Nucl. Phys. A516, 429 (1990). 

[14] J. I. Kapusta, Finite- Temperature Field Theory (Cambridge University Press, 
Cambridge, England, 1989). 

[15] M. LeBellac, Thermal Field Theory (Cambridge University Press, Cambridge, Eng- 
land, 1996). 

[16] M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, Prog. Theo. Phys. 108, 929 
(2002). 

[17] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields (McGraw-Hill Book 
Co., New York, 1965). 

[18] J. Wess and J. Bagger, Super symmetry and Supergravity (Princeton University 
Press, Princeton, 1992). 



17 



